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We use the Relativistic Mean Field (RMF) method to calculate properties of neutron 
rich, usually deformed nuclei, important for equation of state calculations and which 



have significant abundance in the early stages of stellar collapse. We compare the results 
of our microscopic calculations with existing cold nuclear equations of state based on 



macroscopic liquid drop model and the FRLDM model. 



1. Introduction 
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a: 

Spectroscopy of neutrinos produced by electron capture on neutron-rich f-p shell nuclei 
having significant abundance in the cores of nearby presupernova and collapsing stars, 
and emitted before neutrino trapping sets in (at core density of ~ 10 12 gm/cm 3 ), can yield 
useful information on the physical conditions and nuclear composition of the core because 
they stream freely through the overlying stellar matter without any further interactions^ . 

The neutrino spectrum depends, among other factors, on nuclear abundances. The 
post-silicon burning (p > 10 9 gm/cm and T > 6 x 10 9 deg. K) core composition con- 
sists of neutron rich f-p shell nuclei whose abundance is controlled by nuclear statistical 
equilibrum. The nuclear abundances thus depend on the temperature dependant nuclear 
binding energies, and through the neutron and proton fractions X n and X p , on the neutron 
and proton chemical potentials /i n and fi p . The mass of the homologously collapsing core 
and the strength of the hydrodynamic shock after stellar core bounce are determined by 
the lepton fraction and consequently by the electron capture rates®. Since the e~-capture 
threshold too, depends on the nuclear binding energies it is neccessary to calculate the 
chemical potentials and binding energies accurately, as functions of core temperature and 
density. 

Because the temperature and the drip neutron fraction are relatively low, the cold nu- 
clear equation of state (EOS) used should reproduce well the laboratory values of nuclear 
binding energies and chemical potentials for the neutron rich nuclei of interest. Since 
it must also take into account the nuclear shell and pairing effects which persist upto 
temperatures ~ 0.5 to 1 MeV, these quantum effects have to be calculated using mi- 
croscopic mean-field Hartree Fock or RMF methods. Thus, EOS based on the classical 



*NAS-NRC Senior Research Associate; on leave of absence from Tata Institute of Fundamental Research 



2 



liquid-drop modelscl~0, can only hold accurately beyond ptrapping — 10 12 gm/cm 3 and T 
~ 1 MeV where nuclear shell and pairing effects are washed out. For physical conditions 
before neutrino trapping sets in, the cold EOS should go smoothly into the high tem- 
perature macroscopic liquid drop model based EOS0. Since analytical forms of EOS are 
less computationally expensive than detailed microscopic calculations when incorporated 
into hydrodynamic codes to simulate stellar collapse and explosion, we look for analytical 
models which would compare well with the microscopic RMF calculations. 



2. Existing Equations of State 

In general, all EOS developed so far approximate the ensemble of heavy nuclei and drip 
nucleons by a lattice consisting of a single heavy spherical nucleus (A, Z) (which is the 
most bound system for a given core configuration) immersed in a sea of drip neutrons of 
density n n and electrons of density n e . 

In the cold, incompressible liquid drop formalism, the nuclear matter energy of a single 
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nucleus in a lattice, Wn, is given 

by §,§ ; 

W N (x,p N , V N ,u) = W bulk + 290x 2 (l - x) 2 Al + (3x 2 p 2 N V^(l - ~u* + ~u) (1) 
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where x is the proton fraction , p^ the nuclear density, Vjv the nuclear volume and u 
the fraction of total volume occupied by nuclei (u = p/pn)- This model accounts for the 
effects of drip neutrons and the nuclear lattice, but not for the nuclear deformation or 
shell and pairing effects. Minimising the energy of the system with respect to the nucleons 
bound in nuclei at fixed n^A, ri^Z, n n and u^Vn gives the mass number of the (most 
bound) 'mean' nucleus A. 

3. Results 

We used the microscopic RMF approach with the Lagrangian set with non-linear self 
interactions for the cr-meson to calculate nuclear properties for isolated nuclei at zero 
temperatures, since it has been known to reproduce the ground state properties of /in- 
stable nuclei with sufficient accuracy. The parameters that enter into the Lagrangian 
include the nucleon mass M B , the masses of the a, u and p mesons (m CT , m w , m p ) and 
the coupling constants g a g w and g p . These are self consistently determined by variational 
calculation. The mesonic masses and coupling constants are treated as fixed parameters, 
their values having been derived by fitting to ground state properties of a few select 
spherical nuclei. The pairing is dealt within the BCS approximation and the initial values 
of the pairgaps were obtained from the odd-even mass differences, but it was found that 
to reproduce experimental binding energies, it was neccessary to decrease pairgaps as the 
nuclear assymetry increased. Fig.(l) compares the experimental binding energies with the 
RMF values for a range of nuclei in interest to this stage of collapse. Nuclear deformations 
are shown in Fig. (2) for spherical Ni and deformed Zn isotopes. The neutron chemical 
potential p n is obtained from the incompressible liquid drop model asfl: 

a n = -16 + 125(0.5 - x) - 125(0.5 - x) 2 - 290x 2 (l - x) 2 A^ ^ - 7x ] (2) 

2(1 — x) 



Figure 1. Binding Energies vs. Neutron no. N for Mn, Fe, Co, Ni and Cu isotopes 

Figure 2. Quadrupole(/3 2 ) deformation vs. Neutron no. N for Ni and Zn isotopes 

Figure 3. Comparision of analytical models with RMF calculations for Ni isotopes 

Figure 4. Comparision of analytical models with RMF calculations for Zn isotopes 



4 



The first 3 terms in this expression come from the nuclear bulk energy, the last from 
coulomb and surface contributions to the nuclear energy. 

Since the nuclear matter is being treated as an ensemble of isolated cold nuclei with a 
low drip neutron fraction, we can also derive a nuclear equation of state using the Finite 
Range Liquid Drop Modell (FRLDM). This gives the following expression for the volume 
contribution to fi n : 

Vn\voi = -16.00 + 123.04(0.5 - xf + 246.08x(0.5 - x) (3) 

and an expression similar to that in eq.(Q) for the finite size and coulomb effects. Note 
that the bilinear term in eq.(|3|) reduces to the linear term in eq.(|2|) when nuclei in neigh- 
bourhood of x ~ 0.5 are considered. The results of the FRLDM model are compared with 
the RMF and values from eq.(Q) in Fig. (3) for the spherical Ni systems and in Fig. (4) for 
the deformed Zn systems. We find that the FRLDM model reproduces better the RMF 
values than eq.(0) even in the case of spherical Ni nuclei. 

3.1. Extension of the EOS to finite temperatures. 

At finite temperatures, the free energy F of the system is modified to@l 

F = W N -~ T 2 (4) 
A m 

where a is the nuclear level density parameter, and m* is the effective nucleon mass, 
which in general is a function of the nuclear temperature. We are using the RMF code 
to calculate the level density of the last filled orbitals to extract the coefficient of the 
temperature dependent correction. 

This work was carried out under 8th Five Year Plan Project 8P-45 at T.I.F.R. 
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